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The problem of hnding a marked node in a graph can be solved by the spatial search algorithm 
based on continuous-time quantum walks (CTQW). However, this algorithm is known to run in 
optimal time only for a handful of graphs. In this work, we prove that for Erdos-Renyi random 
graphs, i.e. graphs of n vertices where each edge exists with probability p, search by CTQW is almost 
surely optimal as long as p > log^^^(n)/n. Consequently, we show that quantum spatial search is 
in fact optimal for almost all graphs, meaning that the fraction of graphs of n vertices for which 
this optimality holds tends to one in the asymptotic limit. We obtain this result by proving that 
search is optimal on graphs where the ratio between the second largest and the largest eigenvalue is 
bounded by a constant smaller than 1. Finally, we show that we can extend our results on search to 
establish high fidelity quantum communication between two arbitrary nodes of a random network of 
interacting qubits, namely to perform quantum state transfer, as well as entanglement generation. 
Our work shows that quantum information tasks typically designed for structured systems retain 
performance in very disordered structures. 

PACS numbers: 03.67.Ac, 03.67.Lx, 03.67.Hk 


Quantum walks provide a natural framework for tackling 
the spatial search problem of finding a marked node in a 
graph of n vertices. In the original work of Childs and 
Goldstone it was shown that continuous-time quan¬ 
tum walks can search on complete graphs, hypercubes 
and lattices of dimension larger than four in 0{y/n) time, 
which is optimal. More recently, new instances of graphs 
have been found where spatial search works optimally. 
These examples show that global symmetry, regularity 
and high connectivity are not necessary for the optimal¬ 
ity of the algorithm US]. However, it is not known how 
general is the class of graphs for which spatial search by 
quantum walk is optimal. 

Here we address the following question: If one picks at 
random a graph from the set of all graphs of n nodes, 
can one find a marked node in optimal time using quan¬ 
tum walks? We show that the answer is almost surely 
yes. Moreover, we adapt the spatial search algorithm 
to protocols, for state transfer and entanglement genera¬ 
tion between arbitrary nodes of a network of interacting 
qubits, that work with high fidelity for almost all graphs, 
for large n [nodes and vertices are used interchangeably 
throughout the paper). Thus, besides showing that spa¬ 
tial search by quantum walk is optimal in a very general 
scenario, we also show that other important quantum in¬ 
formation tasks, typically designed for ordered systems, 
can be accomplished efficiently in very disordered struc¬ 
tures. 

We obtain our results by studying the spatial search 
problem in Erdos-Renyi random graphs, i.e. graphs of 
n vertices where an edge between any two vertices ex¬ 
ists with probability p independently of all other edges, 
typically denoted as G[n,p) [S] [5]. Note that our ap¬ 
proach is different from the quantum random networks 


of non-interacting qubits defined in [7], where two nodes 
are connected if they share a maximally entangled state, 
having in view long-distance quantum communication. 
Also, in Refs. BM, the authors compare the dynamics 
of classical and quantum walks on Erdos-Renyi graphs 
and other complex networks, although with a different 
perspective from our work. 

In our work, we show that search is optimal on G[n,p) 
with probability that tends to one as n tends to infinity, 
as long as p > log^^^(n)/n. It can be demonstrated that 
when p = 1/2, G(n, 1/2) is a graph picked at random 
from the set of all graphs of n nodes in an unbiased way, 
i.e. each graph is picked with equal probability. This al¬ 
lows us to conclude that spatial search by quantum walk 
is optimal for almost all graphs from this set. To obtain 
this result, we prove a sufficient condition regarding the 
adjacency matrix of graphs where search is optimal: the 
eigenstate corresponding to its largest eigenvalue must be 
sufficiently delocalized and the ratio between the second 
largest and the largest eigenvalues must be bounded by 
a constant smaller than 1. 

This general result also allows us to prove that search 
is optimal for graphs sampled uniformly from the set of 
all regular graphs, also known as random regular graphs. 
Thus, this leads us to conclude that spatial search by 
quantum walk is optimal for almost all regular graphs. 

A sufficient condition for optimal quantum 
search - Let G be a graph with a set of vertices 
V = {!,..., n}. We consider the Hilbert space spanned 
by the localized quantum states at the vertices of the 
graph PL = span{|I),..., |n)}, and the following search 
Hamiltonian 


Hg = - |w) (w| - uAg, 


( 1 ) 


2 


where |ui) corresponds to the solution of the search prob¬ 
lem, 7 is a real number and Ac is the adjacency matrix 
of a graph G [T] . We say that quantum search by contin¬ 
uous time quantum walk is optimal on a graph G if there 
is an initial state |V’o)) irrespective of w, and a value of 
7 such that after a time T = 0{^/n) [10], the probability 
of finding the solution upon a measurement in the ver¬ 
tex basis is | {w\e~^^^*\tljo) p = 0(1). The initial state 
|'0o) is usually chosen to be the equal superposition of 
all vertices, i.e. the state |s) = K) /v^j since it is 

not biased towards any vertex of the graph. We start 
by proving the following general lemma regarding the 
spectral properties of Aq and the optimality of search: 

Lemma 1 Let Hi be a Hamiltonian with eigenvalues 
Ai > A 2 > ... > Afe (satisfying Ai = 1 and |Ai| < c < 1 
for all i > 1) and eigenvectors |?;i) = |s), \v 2 ) ,■■■, \vk) 
and let H 2 = |w) (w| with |(i(;|s)| = e. For an ap¬ 
propriate choice of r = 0(1), applying the Hamiltonian 
(1 -I- r)Hi H 2 to the starting state |ni) = |s) for time 
0(l/e) results in a state \f) with Kw |/) p > ~ o(l). 

Proof: See Section I in Supplemental Material. 

Thus, if Af > A^ > ... > Aj^ are the eigenvalues of 
the adjacency matrix Aq, we choose 7 = 1 /Aj^ and con¬ 
sequently, Hi = jAq. If |s) is an eigenstate of Ac 
corresponding to its largest eigenvalue A)^, and since 
Kw |s) I = Ify/n, we have that search is OTtimal as long 
as A^/Aj^ < c < 1, following Lemma We will see 
that Erdos-Renyi graphs and random regular graphs ful¬ 
fil this property, leading to the conclusion that search 
is optimal for almost all graphs and also for almost all 
regular graphs (the latter is discussed in Section II of 
Supplemental Material). 

In fact, Lemma [^implies that for any regular graph hav¬ 
ing a constant normalized algebraic connectivity, quan¬ 
tum search is optimal HI]. This is in contrast to Ref. gj 
where two examples of regular graphs |I 2 ] with low nor¬ 
malized algebraic connectivity are given, such that quan¬ 
tum search is optimal on one and non-optimal on the 
other. This result showed that normalized algebraic con¬ 
nectivity is not a necessary condition for fast quantum 
search: when connectivity is low, search can be fast or 
slow depending on the graph. On the other hand. Lemma 
proves that high connectivity is indeed a sufficient con¬ 
dition. 

Quantum search on Erdos-Renyi random graphs 

- Let us consider a graph G{n) with a set of ver¬ 
tices V = {I,...,n}. We restrict ourselves to simple 
graphs, i.e. graphs which do not contain self-loops or mul¬ 
tiple edges connecting the same pair of vertices. The 
maximum number of edges that a simple graph G{n) 
can have is iV = ( 2 ). Thus, there are (^) graphs of 
M edges and the total number of (labelled) graphs is 

J2m=o (m) “ [131 ■ consider the random 

graph model G{n,p), a graph with n vertices where we 
have an edge between any two vertices with probability 
p, independently of all the other edges ISllSlIIlj. In this 
model, a graph Go with M edges appears with proba¬ 


bility P{G(n,p) = Go} = p^(I — p)^~^. In partic¬ 
ular, if we consider the case p = 1/2, each of the 2^ 
graphs appears with equal probability P = 2~^. In 
their seminal papers, Erdos and Renyi introduced this 
model of random graphs and studied the probability of 
a random graph to possess a certain property Q mi. 
They studied properties like connectedness of the graph, 
the probability that a certain subgraph is present, etc. 
They introduced the terminology stating that almost all 
graphs have a property Q if the probability that a random 
graph G(n,p) has Q goes to I as n —>■ 00 . Equivalently, it 
can be stated that G{n,p) almost surely has property Q. 
Interestingly, certain properties of random graphs arise 
suddenly for a certain critical probability p = Pc, where 
this probability depends typically on n. More precisely, 
if p(ji) grows faster than pdn), the probability that the 
random graph has property Q goes to I in the asymp¬ 
totic limit, whereas if it grows slower than Pc{n) it goes 
to 0. For example above the percolation threshold, i.e. 
when p > log(n)/n the graph is almost surely connected, 
whereas if p < log(n)/n the graph has almost surely iso¬ 
lated nodes. 

In this work, we are interested in the threshold value 
of p for which quantum search becomes optimal, i.e. a 
marked vertex from the graph can be found in 0{ySi) 
time. We consider the search Hamiltonian in Eq. 0 
for Erdos-Renyi random graphs Hq(^p^ = — |w) (w| — 
')Ac(n,p)- In order to apply Lemma [T| we need to know 
the largest eigenvalue of AG{n,p)i which we denote as Aj^, 
its corresponding eigenstate |ui) and the second largest 
eigenvalue of ^G(ra,p) denoted as A^. It was shown in 
Ref. [T3] that the highest eigenvalue, Xf is a random vari¬ 
able whose probability distribution converges to a Gaus¬ 
sian distribution with mean np and standard deviation 

as n —>■ 00 . The corresponding eigenstate, 
|wi) tends almost surely to |s) = I/'\/IISr=i I*)- ^ 

more detailed analysis of the convergence of |?;i) to |s), 
refer to Lemma 2 in Section III of Supplemental Mate¬ 
rial. It is also possible to obtain an upper bound on the 
second highest eigenvalue, A^ from the results of Ref. m 
which applies to random symmetric matrices. In fact in 
Ref. [T^, a tighter bound on A^ is provided as n — 00 , 
given by 

X 2 = 2y/np + 0{{npY^^ log(n)) (2) 

We see that as long as p > log^'^^(n)/n, the ratio A^/A(^ 
is bounded by a constant. However, as can be seen in 
Section HI of Supplemental Material, in order to en¬ 
sure that |ui) converges to |s), almost surely, we choose 
the critical value of probability for search to be optimal 
as p > log^^^(n)/n. In fact, in the asymptotic limit, 
A^/Af —>■ 0, and the eigenstates corresponding to the 
two lowest eigenvalues of HG(n,p) are 

IX \ _ 1 ^) ± /ox 

|A±)-(3) 

where |su)) is the equal superposition of all the vertices 

other than the solution state jic). The probability of 




p=0.1 
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success is 

Pw{t) = I (w| exp{-iHa(n,p)t)\s) = sin^ (;^) ’ 

To confirm these theoretical predictions we plot, on the 
left side of Fig.[T]a)-c), the approximate probability Pwit) 
from Eq. ffl (in red) and the exact solution calculated nu¬ 
merically (in blue) for n = 1000 and p = 0.1,0.01,0.002. 
On the right side, we plot the spectrum of the respective 
Hamiltonians. We observe, as expected, that the larger 
the gap between the two lowest eigenvalues and the bulk 
of the spectrum, the better is the approximation given 
by Eq. Q for the probability of success of search. As 
this gap disappears, close to the percolation threshold, 
the eigenstates corresponding to the two lowest eigenval¬ 
ues do not follow Eq. (§ and will mix randomly with 
the subspace orthogonal to \w) and Ists)- At this point, 
since we are close to the percolation threshold, the graph 
is expected to have some isolated components and the 
algorithm breaks (see Fig. He)). 

So far we have made the choice 7 = 1 /A(^, and assumed 
that we know the value of the random variable In 
fact, its standard deviation is small enongh so that it is 
snfhcient to know its mean, which is equal to np, i.e we 
can choose 7=1/ {np) , in order to prove that search 
is optimal almost surely. We prove this in Section IV 
of Supplemental Material, using tools of degenerate per¬ 
turbation theory. These tools are also useful to design 
protocols for performing optimal state transfer and en¬ 
tanglement generation in Erdos-Renyi graphs, as will be 
explained subsequently. 

State transfer with high fidelity - Quantum state 
transfer in spin chains Hi! and spin networks [ 18 ] has 
been proposed as a way to establish short-range quan¬ 
tum channels. The problem of what structures lead to 
high fidelity state transfer has been of wide interest (TS)- 
[20] . Here we show that it is possible to transfer, with 
low control and high hdelity, a quantnm state between 
two arbitrary non-adjacent nodes of a random network 
(namely, an Erdos-Renyi random graph). The Hamil¬ 
tonian of a network of coupled spins, with an XX type 
interaction, conserves the number of excitations and so, 
in the single excitation subspace, the Hamiltonian is that 
of a single particle quantum walk on the same network. 
The graph G{n,p) can be perceived as a communication 
network where each node represents a party that trans¬ 
fers information to any of the other nodes. We assume 
that each party has access to a qubit and can control 
the local energy of the corresponding node. In order to 
transfer a state from node i to j, with fidelity that tends 
to 1 in the asymptotic limit, the strategy is the following: 
all qubits are initially in state |0), which is an eigenstate 
of the network; the sender (corresponding to node i) and 
the receiver (corresponding to node j) can tune the re¬ 
spective site energies of |*) and |j) to —1, thereby making 
|z), |j) and |s) approximately degenerate. Finally, in or¬ 
der to transfer a qubit from i, the sender performs a local 
operation on her qubit to prepare |V’) = a | 0 ) -I- /? | 1 ). As 
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FIG. 1 : Left side: probability of observing the solution cal¬ 
culated numerically (blue curve) compared to the prediction 
from Eq. (red curve), obtained in the limit n —>■ 00, using 
degenerate perturbation theory. We hx the number of vertices 
n = 1000 and p = 0.1, 0.01 and 0.002 in a), b) and c), respec¬ 
tively. Right side: Spectrum of the search Hamiltonian for 
instances of random graphs that provide the dynamics repre¬ 
sented on the left side. In red, the two lowest eigenvalues are 
shown in a) and b), which are clearly isolated from the rest 
of the spectrum shown in blue. In c) this does not happen 
since p is close to 1/n, which is the percolation threshold and 
thus the semicircle law is not valid. We see that, the larger 
the gap between the two lowest eigenvalues A± and the rest 
of the spectrum, the better is the prediction from Eq. for 
the probability of success. When the two lowest eigenvalues 
are not isolated, the probability of observing the solution is 
low and the algorithm does not provide speed-up with respect 
to classical search. 


long as p > log^^^ (n) /n, the approximate dynamics of a 
quantum walk starting at |z) is obtained by diagonalizing 
the Hamiltonian 


^G(n.p) = - K) (*l - \j) 01 - Oij) {Sfj\ - l^G{n,p) (5) 


projected onto the approximately degenerate subspace 
spanned by {|*), jsjj), |j)} which is given by 


^G{n,p) 


-1 -1/v^ 0 

-ll^/n -1 _ -1/v^ , 

0 - 1 / 0 ^ -1 


( 6 ) 
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FIG. 2 : Quantum state transfer in Erdos-Renyi random 
graph G( 100 , 0 . 2 ): using our protocol, the fidelity achieved 
for this network is 80 %. 


with |s-) = Nn-2 and |s-) « |s-) «Js), 

where we assume that i and j are non-adjacent vertices. 
Thus, the dynamics is approximately the same as that 
of end-to-end state transfer in a chain with three spins, 
where perfect state transfer is possible [5D] and the com¬ 
ponent of the wave function at the receiver is approxi¬ 
mately I {j\U{t)\i) p = Hence, after time 

T = 'iT^Jnl2, the receiver gets IV') with fidelity 1, in the 
limit n —>■ oo (see Fig. |^for an example with finite n). 
The receiver can preserve this state for future use by 
tuning the energy of node j, locally, to a value that is 
off-resonant with the rest of the network [5T]. We con¬ 
clude that high fidelity quantum state transfer can be 
achieved in almost all networks. 

Creating Bell pairs in a random network - In quan¬ 
tum communication networks, entanglement is an useful 
resource that can be used for various tasks such as tele¬ 
portation, superdense coding, cryptographic protocols, 
etc [22]. Here, we present a protocol to entangle arbi¬ 
trary nodes in a random network based on the search 
Hamiltonian. Imagine that Charlie at node |ui) wants 
to entangle the qubits of Alice at node |o) and of Bob 
at node \b). We assume that none of the nodes |w), |a) 
and \b) are adjacent to each other. As before, 7 is cho¬ 
sen to be l/{np). In this case, the protocol is as fol¬ 
lows: i) Alice, Bob and Charlie tune their respective site 
energies to —1, ii) Charlie tunes his nearest neighbour 
couplings to y/2ldc, where dc is the degree of the node 
corresponding to Charlie, while the other couplings in 
the graph are 7 = 1/np. This ensures that the Hamilto¬ 
nian, projected onto the approximately degenerate sub¬ 
space spanned by \w), |s^) = 'Ek^a,b.w |fc) and 

|sab) = (|a) + \b))/V^, is equal to 


^G{n,p) 


-1 0 
-^/2jn -1 -y/Vn , 

0 -V^ -1 


(7) 


in the asymptotic limit [23] . Thus, after time T = 
7r-yn/2, Alice and Bob share the state [Sab) = (|a) -f 


\b))/y/2, which is a Bell state. Subsequently, other Bell 
states may be obtained by local operations. Furthermore, 
Alice and Bob can preserve their Bell state by tuning the 
local energies of their qubits to a value that is off-resonant 
with the other eigenvalues of the network. 

Discussion - We have shown that searching for a 
marked node in a graph using continuous-time quantum 
walks works optimally for almost all graphs. This means 
that, in terms of the structures on which it performs op¬ 
timally, this approach to quantum spatial search is much 
more general than what has been shown before. Our re¬ 
sult was obtained by proving that the algorithm is almost 
surely optimal for Erdos-Renyi random graphs G{n,p), 
as long as p> log^^^(n)/n. 

As pointed out in Ref. [T] , the analog version of Grover’s 
algorithm of Ref. m can be seen as a quantum walk on 
the complete graph. Furthermore, Erdos-Renyi random 
graph G(n,p) can be obtained from the complete graph 
by randomly deleting edges with probability 1 —p. Thus, 
our result can also be interpreted as showing an inherent 
robustness of the analog version of Grover’s algorithm to 
edge loss. This implies that there is a large family of 
random Hamiltonians that can be employed to achieve 
optimal quantum search. Hence, our work paves the way 
to understanding how this randomness would translate 
to the circuit model of quantum search and whether this 
implies an inherent robustness of the (standard) Grover’s 
algorithm. 

Finally, we have shown that one can adapt the spatial 
search algorithm to design protocols for quantum state 
transfer and for entanglement generation between arbi¬ 
trary nodes of a random network of interacting qubits. 
Our results show that quantum information tasks typi¬ 
cally designed for structured systems retain performance 
in very disordered structures. These results could lead to 
further investigation on what kind of random structures 
appear naturally in physical systems (for example those 
appearing in Refs. [211 [21]) and whether they would offer 
a sufficient spectral gap to perform efficient and robust 
quantum information tasks. It would also be interest¬ 
ing to explore whether non-trivial quantum information 
tasks can be performed on other models of random net¬ 
works such as scale-free networks [26j . 
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SUPPLEMENTAL MATERIAL 

I. SUFEICIENT CONDITION FOR OPTIMAL 
QUANTUM SEARCH: PROOF OF LEMMA 1 


Proof: Let us express 

\w) = ai |ui) + 02 \V 2 ) - 


']) in the basis 


■ flfe \Vk) ■ 


( 8 ) 


We rescale Hi by (1 + r)Hi — rl, where I is the identity 
matrix. With this replacement, |ui) remains an eigen¬ 
vector of Hi with an eigenvalue 1. The other eigenvalues 
change to A' = (I -I- r)Ai — r. Now, the expression 


E 

i=2 


1-A,’ 


(9) 


after rescaling Hi as mentioned before becomes 

k n k 9 




(1 -I- r)(l — Ai) 


( 10 ) 


As I Ail < c, we have 


E 


^E 


a? 


i^2 (1 +^(1+0(1 -A) 


, and. 


E 

2=2 


(1 -t- r)(l — \i) 


^E 

2=2 


a? 


(l-f r)(l-c) 


( 11 ) 

( 12 ) 


So we choose an appropriate r S [—yfy:, yiy;], such 
that 




(13) 


i=2 i=2 

If |Ai| < c, then A' > 0 for r = —yfv and A' < 0 for 


r = In the first case, the left hand side of (13) 

is at least the right hand side. In the second case, tile 
left hand side is at most the right hand side. After the 
replacement of Hi by (1 -|-r)7Li — r/, the new eigenvalues 
A 2 ,..., AJj are in the interval [—i^]- 
After we replace Hi + H 2 by (1 -|- r)Hi + H 2 — rl, we 
can omit the —rl term (since it only affects the phase of 
the state). To simplify the notation, we now refer to the 
Hamiltonian (1 -|- r)Hi as Hi and to new eigenvalues A' 
as Ai. 

Let 

|w) = bi |ui) -I- &2 1 ^ 2 ) + ... + bk \vk). (14) 

We write out the conditions for |u) to be an eigenvector 
of Hi + H 2 with an eigenvalue A. We have 

H |u) = Hi |u) -f H 2 |u) = ^(6iAi -f Oiy) |uj) (15) 


where 7 = {w\v). Since we also have 

H\v) = A|u) = ^A 6 , |ui) 


(16) 


we get that A&i = Xibi + yui which is equivalent to bi = 
. Substituting this into 7 = (w |u) = J2i O'ibi gives 


A-Ai 

that 


E 


A-A 


= 1 . 


(17) 


This is the condition for the eigenvalues A. In each of 
intervals [Ai, Ai_i] ioi i = 2,k, the left hand side is 
strictly decreasing from -l-oo to —00 and in the interval 
[Ai,-|-oo), the left hand side is strictly decreasing from 
-foo to 0. Therefore, each of these intervals contains one 
eigenvalue. 

We are interested in the two eigenvalues A that are in 
[A 2 ,Ai] and [Ai,-|-(X)). (We denote these eigenvalues by 
A_ and A+ and the corresponding eigenvectors by |z;_) 
and |u+).) We express these eigenvalues as A = 1 -f <5 
(where S is positive for A+ and negative for A_). By 
Taylor expansion, if 6 is small, we have 


E 


(l + A-A 


_ _ O'i _ _ 0 

- X- “ ^ 1 - A,; ^ - 


( 1 -A)^ 


Thus, the condition for eigenvalues becomes 


5 + 0 (, 52 ). 

(18) 




i=2 


1-A,; 


i=2 


(1 - A) 


Since the second term on the left hand side is 1, this is 
equivalent to 


^ = V_ - _ 5 + 0(6'^) 

A Z^(1-A)2 1 ^ 


i=2 

which is satished for 

5 « ± 


Oi 


( 20 ) 


( 21 ) 


E k a 

713 


(1-A.) = 


Since oi = e and the denominator is of the order 0(1), 
the right hand side is 0 (e). 

We now consider the overlap (s|z;+). We assume that 
|u+) is normalized so that ||z;+|| = 1. This is equivalent 
to = 1 which, in turn, is equivalent to 


E 


7 

A — Ai 


= 1 . 


( 22 ) 


We can rewrite this as 
1 


7 




(A-A)2- 


(23) 
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We now estimate the expression under the square root. 
We have 


1 


7 \ 5^ ^ ^ 


(24) 


with the first step following by approximating A — = 

(1 + 5) — Ai « 1 — Ai for i > 1 and the second step follows 
from ( 21 ). 

This means that 7 « —. Therefore, 

v 2 ai 


(s \v+) = 


701 

A ^-1 


7«i _ 1 

S y /2 


and 


u_) = 


701 

A_ - 1 


7^1 

(5 


1 

'Ti' 


(25) 


(26) 


Thus, |s) can be approximated by ^(|z;+) —|z)_)). Evolv¬ 
ing the Hamiltonian Hi + H 2 for time ^ = 0(^) trans¬ 
forms |s) to 




(27) 


We have 


{w\f) --^{w\v+) +-^{w\v-). (28) 

We now consider {w |r;+) = 7 . By combining the first 
part of (24) with (21), we obtain that 


2E 




(29) 


Because of (13) and Xi < this is at most 


2 S ^<.2 




2c — 
1+c 


1 - ^ 
1+c 


2 ( 1 -He) 
1 — c 


(30) 


Therefore, 7 = (i(;|i;+) > Similarly, (w |w_) > 

2 ( 1 +c) ■ Together with (28), this means that, up to the 
approximations that we inade. 


'!/)> 


1 — c 
1 -I- c’ 


(31) 


□ 


II. QUANTUM SEARCH ON RANDOM 
REGULAR GRAPHS 

A family of random graphs whose adjacency matrix has 
an 0 ( 1 ) gap between the largest and second largest eigen¬ 
values are the d-random regular graphs, a random graph 


sampled uniformly from the set of all regular graphs of 
degree d. For these graphs, the largest eigenvalue is d, 
with the corresponding eigenvector, |s). Also in Ref. m 
it has been proven that the second largest eigenvalue is 
0(d^/'*) for d > 3, with high probability. This way, we 
choose 7 = 1 /d and since A^/A(^ = 0 (d“^/^) < 1 , it fol¬ 
lows from Lemma[^that quantum search is optimal. It is 
interesting to note that for lattices, the lowest dimension 
for which search is possible in 0{y/n) time is dimension 
five [ 1 ], which is a specihe instance of a regular graph of 
degree ten. However, for random regular graphs search 
is optimal for degree three and larger. 


HI. GONVERGENCE OF THE EIGENSTATE 
CORRESPONDING TO THE MAXIMUM 
EIGENVALUE OF AN ERDOS-RENYI RANDOM 
GRAPH 

Lemma 2 Let A be the adjacency matrix of the Erdos- 
Renyi random graph G{n,p) with vertices 1,2, ...,n. 
Let 7 A represent the adjacency matrix of G{n,p) with 
each entry rescaled by j = 1/np. Also let |s) = 
(l/-\/n) |i) be the equal superposition of all nodes 

such that |s) = a |z)i) -I- /3 jui)'*", where |r;i) is the eigen¬ 
vector corresponding to the highest eigenvalue, Ai of^A. 
Then, a > 1 — o(l) almost surely for p > _ 

Proof: First we observe that from Ref. (TS] that the 
largest eigenvalue follows a Gaussian distribution with 

mean 1 and standard deviation i.e, Ai 

n y p ‘ ‘ 

A/'(l, So if (5 = \l^fn, for p > log^^^n/n, one 

can show that 


Pr[Xi > 1 — (5] = 1 — 


log^/'‘ n 


V2 


where erfc[x] = ( 2 /y^) e As erfc[x] 

n —>■ 00 , we have 


(32) 
0 as 


Ai ^1 — S, 


(33) 


almost surely. To prove this explicitly, we use the 

_^2 

bound, erfclxl < - 7 , , for x > 0. In our 

’ ^ J (a;+yT2+4/ll) ’ 

case X = (log^/^ n)l\/2 and so, by using the inequality 
log“ n > alogn, for a > 1 , we can show that 


erfc[x] < o(- 


l3/4 log 


3/4 


:)^ 


Thus, 


Pr[Xi >l-6]>l-0(- 

V 7 


l3/4 log 


3/4, 


(34) 


(35) 


Similarly, one can also obtain an upper bound Ai < l-|-d, 
almost surely. 
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Also, from Ref. nsun], we have, 

|| 7 (A-£’(A))|| < (2 + (np)"^/'‘logn)-^, (36) 


where are eigenvalues of the adjacency matrix Aq |28j . 
In the limit of n —> oo this approaches a continuous func¬ 
tion. For a random graph G{n,p), as long as np —>■ oo, 
the spectral density is given by 


where E{X) denotes the expectation of random variable 
X and II II denotes the standard Euclidean norm. 

Let Xi,j > 2 be the rest of the spectrum of 7 A and let 
|ui) be the corresponding eigenvectors. From Eq. 36 it 
follows that 


||Ai |ui) (ui| \vi) {v^\ - |s) (s| II (37) 

i>2 

< (2-b (np)"^/'‘logn)^= (38) 

Now, 

(Ai |ui) (ui| -b ^ Ai |u*) {vi\- |s) (s| ) |ui) (39) 

i>2 

= Ai|ui)-a|s) (40) 

= (Ai - a^) |ui) - a/3|ui)-^ . (41) 


So, 

||(AiK)(ui|+^AK)(ud-|s)(s|)K)|p (42) 

i>2 

= (Ai-a2)2 +0,2/32. (43) 

From Eq. and the Cauchy-Schwarz inequality, we ob¬ 
tain 


2n 2 ^ (2 + (np) i/Mogn)2 
[Ai — a ) S - 


> Ai — 


2 


np 

(2 + logn) 


y/Ep 


>a^> 1 - 0 ( 1 ), 


(44) 

(45) 

(46) 


1 3/2 

for p > ^ , where the final expression for a follows 

from the fact that a G [ 0 , 1 ] and that Ai > 1 — i5. 


□ 


IV. PROOF OF OPTIMALITY OF SEARCH ON 
ERDOS-RENYI RANDOM GRAPHS USING 
DEGENERATE PERTURBATION THEORY 

Here, we present an intuitive way to prove the optimal¬ 
ity of search for Erdos-Renyi graphs, G(n,p), by using 
degenerate perturbation theory. This proof is instru¬ 
mental in constructing protocols for optimal state trans¬ 
fer and entanglement generation in these random net¬ 
works. 

The spectral density of a graph G is defined as 

1 "■ 

p(A) = -^<5(A-A) (47) 


{ ^Anp{l-p) - A2 
pA) — N 27rnp(l — p) 

lo 


if |A| < 2^np{l-p) 
otherwise 


(48) 

known as the Wigner’s semicircle law. The highest eigen¬ 
value, Ai, of is isolated from the bulk of the spec¬ 

trum and follows a Gaussian distribution with mean np 
and standard deviation ^Jp{l — p), as n —> 00 . From 
Section I, the corresponding eigenstate, |ui) tends almost 
surely to |s) = l/-\/nSr=i I*)- From Ref. [16], we obtain 
that the second highest eigenvalue, as n —>■ 00 , is given 
by 


A 2 = 2y/np + 0{(npY^'^ log(n)). (49) 


There is thus a significant gap in the spectrum between 
the first and second largest eigenvalues, with high prob¬ 
ability. In order to make use of this separation between 
the largest eigenvalue and the rest of the spectrum, it 
will be helpful to write 


X^G{n,p) — E |s') (s'l -b (50) 

where E -G np, and |s') —>■ |s) as n —>■ 00 . The search 
Hamiltonian then becomes 


Hg{u,p) = - |w) {w\ - jpE |s') (s'l - 7 p^g(„,p)- (51) 

We use |s) as the initial state of the quantum algo¬ 
rithm and choose 7 p = l/{np) so that, for large n, |s) 
and |ty) are approximately degenerate. The spectrum of 
7 pAg(^ p) follows the semi-circle law, where the radius of 
the semicircle is given by 


R = Jp2^/np{l-p) = 2 J . (52) 

y rip 

As long as np ^ 00 , the radius R ^ 0. This implies, 
for the whole range in which the semi-circle law is valid, 
that the radius R shrinks as np grows. Also from Section 
HI, we know that as long as p > log^'^^(n)/n, 

1-^ < 7pAi < 1-b (53) 

vn Vn 


almost surely. 

One can show that the algorithm retains its optimality 
as long as this is the case. If 


- 1 - (5 < 7 pAi < -1 -b (5, (54) 

using degenerate perturbation theory, we obtain the 
ground and first excited states of Lfc/n.p) from its di- 
agonalization in the 2 -dimensional subspace spanned by 















Il'u;), |sffi)}, where |su)) = Yh^w N) N'^ “ 1- These 
eigenstates are 

\X+) ^ (55) 
|A_) «-(/r|w) + ^ Is^b)) (56) 

where |sbb) = K) /Vn- 1, n = d/2+ 1 / 52 / 4 + 1/n 

and K = fj? + 1/n. 

Thus, the probability of observing the marked vertex |w) 
is given by 

Pw{t) = I {w\exp{-iHG(n,p)t)\s) p = ^ nS^/4 

(57) 


where jA + 1/n. Thus, as long as 5 G 

[—l/-/n, l/-/n], the running time of the algorithm is 

r = o(v^). 


Thus, the gap between the lowest eigenvalue of 
~lp^G(n p) second lowest eigenvalue of HG(^n,p)j 

is 7 p(A 2 — Ai) = 0(1) as long as p > log^^^(n)/n. The 
error obtained from this approximation is since 

7 pA 2 is the largest eigenvalue of ^pA'q^u p) '"^hich we con¬ 
sider as a perturbation. Note that the error of our 
approximation decreases as np increases. Thus, for a 
fixed n, the higher the value of p, the lower is the error. 

’ Again, for a fixed p, the error diminishes with increase in 
n. 
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